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MULTIPLE  SCATTERING  FROM  NONSPHERICAL  TARGETS 


1.  INTRODUCTION 

In  this  report  the  development  of  a  new  extension  of  Shifrin's  scheme  for  calculating  the  scattering 
by  non-spherical,  non-symmetric  particles  is  documented  and  several  applications  presented.  In 
addition  we  have  developed  a  technique  for  properly  computing  the  scattering  of  linearly  polarized  light 
by  randomly  oriented  infinite  cylinders,  as  well  as  finite  cylinders  and  spheroids  and  we  have  applied 
this  technique  to  several  examples. 

The  new  extension  of  Shifrin's  method  involves  taking  targets  of  arbitrary  shape,  but  whose  local 
cross  section  is  circular,  and  breaking  the  target  into  a  series  of  stacked  disks.  This  particular  approach 
permits  the  computation  of  scattering  from  such  targets  as  helices,  wavy  cylinders,  cylinders  whose 
radii  vary,  and  toroids,  to  cite  a  few.  This  new  technique  has  been  tested  against  an  experiment 
conducted  by  Dr.  Ru  Wang  of  the  University  of  Florida  at  Gainsville  on  scattering  of  microwave 
radiation  by  a  dielectric  helix.  The  comparison  between  the  experiment  and  our  calculated  scattering 
patterns  are  included  in  this  report  and  show  good  agreement. 

We  have  also  computed  asymmetry  factors  for  the  scattering  of  polarized  light  from  randomly 
oriented  infinite  cylinders.  This  problem  required  special  care  due  to  the  peculiar  scattering 
characteristics  of  very  long  cylinders.  The  results  of  this  calculation  are  useful  in  computing  the  average 
amount  of  light  scattered  in  the  forward  direction  from  a  collection  of  randomly  oriented  cylinders  such 
as  a  cloud  of  cylinders  or  randomly  oriented  fibers  as  used  in  insulation  materials. 


2.  ASYMMETRY  FACTORS  FOR  THE  SCATTERING  OF  POLARIZED  LIGHT  FROM 
RANDOMLY  ORIENTED  INFINITE  CYLINDERS 

This  computation  is  applicable  to  the  scattering  of  light  from  fibrous  media  such  as  is  used  in 
construction  and  other  applications.  These  fibers  are  characterized  by  a  large  aspect  ratio  (fiber  length/ 
fiber  diameter) .  We  have  recently  shown L2  that  for  finite  fibers  with  aspect  ratios  as  small  as  100,  the 
theory  for  infinite  cylinders  gives  highly  accurate  results.  The  second  of  the  above  references  is  a 
comparison  between  theory  and  an  experiment. 

In  studying  the  passage  of  radiation  through  fibrous  media,  the  main  optical  parameter  is  the  ratio  of 
the  amount  of  radiative  power  scattered  in  the  forward  direction  to  the  power  scattered  in  all 
directions.  In  the  case  of  spherical  particles,  this  is  called  the  asymmetry  factor  p  ,  and  for  spheres 
larger  than  the  wavelength  of  light  (the  size  parameter  x=2;taA  is  greater  than  unity)  it  is  given  by 


p(x)  =Iq1  l(p,x,m)  dp  /  J_i 1  l(p,x,m)  dp  >  1/2 


(11 


where  p  =  cos  0,  with  0  the  scattering  angle,  and  m  is  the  refractive  index. 

For  infinite  cylinders,  equation  (1 )  must  be  modified  to  take  into  account  the  unique  geometry  for 
scattering  from  a  sphere.  Figure  1  shows  this  geometry,  where  the  incident  light  is  along  the  z  axis, 
and  the  target  axis  is  zt.  The  scattered  direction  is  z',  and  the  direction  of  incident  linear  polarization 


Z(  incident) 


Figure  1.  The  geometry  for  the  scattering  process.  The  light  is  incident  along  the  z  axis  and  polarized 
along  the  x  axis.  The  light  is  then  scattered  along  the  z  axis,  and  2<J>0  is  the  angle  between  the  incident 
and  scattered  directions.  The  cylinder  axis  is  denoted  zc,  and  it  s  orientation  angles  with  respect  to  xyz 

are  0  (tilt  angle)  and  v  (  azimuthal  angle). 


defines  the  x  axis.  The  scattering  angle  is  denoted  20o  ,  and  the  orientation  of  the  target  axis  is 
specified  relative  to  xyz  by  the  polar  and  azmuthal  angles  0,  y  . 

For  the  scattering  of  light  from  an  infinite  cylinder,  the  scattering  pattern  is  a  cone  with  apex  at  the 
origin  and  symmetry  axis  along  z  ,  as  shown  in  figure  2.  In  this  figure,  the  line  OA  is  along  the 
direction  of  incidence.  Any  other  element  of  the  cone  is  a  possible  scattering  direction.  For  example, 
the  direction  OC  is  a  scattered  direction  that  is  perpendicular  to  the  direction  of  incidence.  The  point  O' 
is  at  the  center  of  the  base  of  the  cone,  and  the  direction  O’A  is  considered  a  reference  for  the 
measurement  of  a  different  kind  of  scattering  angle  -  it  is  denoted  0  ,  and  it  marks  the  angular  distance 
from  O' A  to  a  point  where  a  scattered  direction  line  intersects  the  base  of  the  cone.  For  example,  the 
angle  Qq  is  the  scattering  angle  for  the  scattered  direction  OC.  It  is  shown  -^that  the  equation  of 
constraint  for  infinite-cylinder  scattering  is  given  by 

sin  0  sin  0/2  =  sin  0O  (2) 

The  scattering  angle  0  for  scattering  in  the  plane  perpendicular  to  the  scattered  direction  is  then  given 
by  (0o  =  450) 

0^  =  |2  sin'1  [1  A/2  sin  0  j  j.  (3) 

We  can  now  properly  define  the  asymmetry  factor  for  scattering  from  an  infinite  cylinder  : 

p(<p.y) = J-ee (li  +l2)de /  J./(h  +l2)de.  (4) 

In  this  equation  Ij  and  I2  are  the  respective  scattering  intensities  for  the  scattering  of  light  respectively 
parallel  and  perpendicular  to  the  scattering  plane  (  see  reference  2). 

Calculations  based  on  the  above  work  are  to  be  found  in  reference  5.  The  asymmetry  factor  of  equation 
(4)  for  the  size  parameter  x=25  and  an  index  of  refraction  m  =  1.5  is  plotted  in  figure  3  against  the 
orientation  angle  0.  The  results  are  strikingly  similar  for  this  wide  range  of  size  parameters  (0.1  to 
50).  Perhaps  more  interesting  is  the  plot  of  the  asymmetry  factor  against  size  parameter  shown  in 
figure  4  where  an  average  over  all  orientations  is  made.  We  see  in  this  case  that  the  factor  seems  to 
converge  to  the  value  0.813. 

3.  APPLICATION  TO  THE  OPTICAL  PROPERTIES  OF  FIBROUS  INSULATION 

This  application  is  distinguished  from  that  given  in  section  2  by  calculation  of  scattering  intensities 
from  a  cloud  or  tenuous  medium  of  long  cylindrical  fibers.  Thus,  we  need  to  define  the  averaging 
procedure,  and  this  turns  out  to  be  complicated  by  the  equation  of  constraint  for  the  scattering  of  light 
from  an  infinite  cylinder  -  the  equation  (2)  of  section  2. 

One  might  think  (  and  previous  calculations  were  based  on  this  idea)  that  an  average  intensity  for 


M 
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Figure  3.  The  scattering  cone  for  iight  scattering  from  a  cylinder  with  tilt  angle  0  >  45  .  The  critical 
scattering  angles  0A  and  -0A  are  shown. 
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the  scattering  of  light  would  be  given  by 

lave  =(l/4rc2)  Jl  1(0, Y)  Sin  p  dp  dy.  (5) 

where  I  is  the  scattering  intensity,  and  p,y  are  the  orientation  angles  for  the  cylinder  axis  as  defined 
in  section  2. 

The  conical  scattering  pattern  from  an  infinite  cylinder  complicates  the  averaging  procedure  -  an 
equation  of  constraint  is  needed.  This  equation  can  be  derived  noting  that  from  first-order  scattering 
theory  (Rayleigh-Gans),  the  scattering  intensity  from  a  finite  cylinder  of  height  h  contains  the  factor 

sin2(ph/2)/(ph/2)2 

where 

P  =  k0  (ko  -  k')  kt . 

Here,  kQ  ,k’,  kt  are,  respectively,  unit  vectors  in  the  incident,  scattered,  and  cylinder-axis 
directions.  As  the  height  of  the  cylinder  increases  to  infinitity,  this  factor  vanishes  unless  p  =  0,  or 

(ko-k*)kt=0. 


This  is,  in  effect,  the  condition  that  must  be  satisfied  for  the  scattering  from  an  infinite  cylinder  to  be 
along  the  cone  (Naturally,  the  exact  theory  for  the  infinite  cylinder  predict  the  conic  scattering  pattern 
and  would  give  rise  to  the  same  condition) .  Using  the  geometry  given  in  figure  1,  this  implies  that 

(k0  -  k’)  kt  =  {1^  -  [cos(2p0)  k0  +  sin(  2p0)  i0]} 

x[sin(0)  cos  (y)  i0  +  sin(p)  sin  (y)  j0  +  cos(p)  k0  ] 

=  0, 
or 

0  =  tan*1  (tan(o0)/  cos(y) )  (6) 

We  want  to  emphasize  at  this  point  that  equation  (6)  was  derived  by  us  ^  and  used  for  the  purpose  of 
summing  over  orientations  for  the  first  time  in  the  literature.  If  we  note  figure  5,  the  orientations  of  the 
cylinder  axis  consistent  with  a  scattering  angle  of  20o  and  equation  (6)  follows  a  great  circle  path  from 
point  A  to  point  B,  and  for  a  unit  sphere,  the  length  of  this  path  is 
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Figure  5.  The  scattering  geometry  shown  in  the  reference  frame  x,y,z,  where  the  x-z  plane  is  the 
scattering  plane.  The  orientation  angles  of  the  cylinder  axis  are  0  and  y,  and  the  polarization  angle  of  the 
linearly  polarized  wave  Eq  is  a.  The  scattering  angle  between  the  directions  of  incidence  and  scattering 
is  denoted  2<>q.  The  scattering  angle  used  in  the  theory'  for  the  infinite  cylinder  is  measured  in  a  plane 
perpendicular  to  the  cylinder  axis,  and  it  is  denoted  0. 


JaB  ds  =  SoK/2  dy  sin2(o)/sin(0o)  =  \q^2  c ty  /(I  -  cos2b0  sin2  Y) 

=  ?t/2. 

Therefore,  the  average  intensity  for  the  scattering  of  light  from  a  cloud  of  randomly 
oriented  long  cylinders  is 

|ave=  (1/2rt)  Jo2*  l(0(Y).Y)  sm2(0)/ sin(0o)  dy  (7) 

The  intensities  to  be  averaged  are  denoted  lj  and  I2  .  as  already  noted,  and  correspond  to  scattering 
in  and  normal  to  the  scattering  plane.  For  an  infinite  cylinder,  these  intensities  are  given  by 

ll  =  const.  |b0|  +  2  I  bn|  cos(n6)  ]  cos  y* 

-i  [2  I  bn||  sin(  n9)  ]  sin  y*  ,  (8) 

12  -  const.  |a0n  +  21  an\\  cos(n6)  ]  sin  y* 

-i[2Ian)  sin(  n0)  ]  cos  y*  ,  (9) 

where  the  "a"  and  "b”  coefficients  are  derived  in  detail  in  reference  7. 

These  coefficients  are  dependent  on  the  size  parameter  of  the  cylinder,  the  index  of  refraction  (which 
may  be  complex),  and  the  tilt  angle  d.  The  angle  y*=  y  -  a,  where  a  is  the  angle  that  the  incident 
electric  field  makes  with  the  scattering  plane. 

We  have  made  three  applications.  The  first  is  a  comparison  with  the  work  of  McKay  et.  al^  in 
which  the  averaging  is  done  using  equation  (5)  instead  of  equation  (7).  We  calculate  the  intensities 
for  light  of  wavelength  10  microns  incident  on  cylinders  with  the  complex  index  of  refraction  m  = 
1.85  -i  0.033  and  size  parameter  x  =  7.3513.  These  are  the  values  used  in  reference  8.  The  results  for 
1 1  and  I2  are  shown  in  figure  6  for  the  incident  light  polarized  in  and  perpendicular  to  the  scattering 
plane  (a=0  and  a=7t/2,  respectively).  Also  shown  in  figure  6  is  the  intensity  I  j 2  where  the  incident 
light  is  polarized  in  the  scattering  plane  and  the  scattered  light  is  polarized  perpendicular  to  the  scattering 
plane. 

The  salient  results  are: 

(1)  A  strong  forward  scattering  value.  This  is  in  contrast  to  the  results  of  reference  8  where  the 
intensities  go  to  zero  in  the  forward  direction  as  a  result  of  compensating  for  a  singularity  that  occurs 
because  of  an  incorrect  averaging  procedure. 

(2)  Smaller  values  in  the  region  of  70  to  1 50  degrees  of  scattering. 

(3)  A  backscatter  value  approximately  1  to  2  orders  of  magnitude  smaller  than  the  forward  scattering 
value. 

(4)  The  number  of  maxima  roughly  equal  to  the  value  of  the  size  parameter. 


Figure  6.  The  relative  intensities  defined  in  equation  9  versus  the  scattering  angle  2<J)g.  The  top  curves 

are  I,  for  a  =  0  and  U  for  a  =  rt/2.  The  bottom  curse  is  1-,  for  a  =  0  (  or  I,  for  a  =  jt/2  ).  The  size 
parameter  x  =  7.35 1 3  and  m  =  1 .85  -  ().()33i  . 
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Figure  7,  The  relative  intensities  defined  in  equation  9  versus  the  scattering  angle  20q  for  x 
17.173  and  m  =  1.29  -  0.0954  i.  a)  The  results  of  reference  3.  b)  our  corrected  results. 
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Figure  9.  The  relative  intensity  defined  in  equation  9  versus  the  scattering  angle  for  x  =  45.7  and  m  = 
1.31.  The  experimental  results  are  shown  by  black  diamonds  and  our  theoretical  results  are  shown  by 
open  squares. 
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The  second  application  is  a  comparison  with  the  work  of  Stephens^  which  also  used  the  wrong 
averaging  procedure  to  determine  the  effect  of  the  scattering  of  thermal  radiation  from  randomly 
oriented  ice  crystals  (  which  are  long  hexagonal  prisms). 

We  calculate  intensities  for  light  of  wavelength  11  microns  .  index  m  =  1.29  -i  0.0954,  and  size 
parameter  x  =  17.136.  The  results  are  shown  in  figures  7  and  8  (ours  and  Stephens,  respectively). 
Our  results  show  marked  differences  especially  at  the  larger  scattering  angles,  where  our  results 
continue  to  show  fluctuations  over  the  full  angular  range. 

Finally,  we  make  a  comparison  of  our  averaging  procedure  with  the  experimental  results  of 
Huffman  et  al.^  In  this  experiment,  visible  light  of  wavelength  0.55  microns  is  incident  on  ice  crystals 
whose  size  parameter  is  x  =  45.7  and  m  =  1.31.  The  results  are  shown  in  Figure  9  where  the 
experimental  points  are  shown  as  black  diamonds  and  the  theory  as  open  squares.  We  see  general 
agreement  over  the  whole  range  of  the  experiment  (for  10  to  150  degrees) .  In  particular,  we  note  that 
agreement  is  achieved  at  the  largest  scattering  angle  of  150  degrees. 

4.  THE  SCATTERING  OF  LINEARLY  POLARIZED  LIGHT  FROM  RANDOMLY  ORIENTED 
SHORT  CYLINDERS  AND  SPHEROIDS 

In  contrast  to  the  first  two  sections,  we  now  study  the  scattering  of  light  from  finite  cylinders  and 
spheroids  of  the  same  material,  equal  volume,  and  the  same  aspect  ratio  (length/diameter).  These 
aspect  ratios  range  from  10  down  to  0.1.  The  purpose  of  this  study  is  to  determine  if  randomizing 
the  target  orientation  effectively  negates  the  differences  in  shape  for  these  targets.  The  iterative 
approach  of  Shifrin  and  Acquista*  ^ is  used. 

The  theory  for  the  scattering  of  light  from  a  finite  cylinder  arbitrarily  oriented  relative  to  the 
incident  light  is  given  in  reference  1  and  this  includes  the  derivation  of  the  second-order  iteration  that 
effectively  includes  the  self-interaction  of  the  various  parts  of  the  target.  An  application  to  finite 
cylinders  and  spheroids  of  equal  volume  and  aspect  ratio  was  done  in  reference  13.  In  reference  13,  the 
second-order  effect  is  evaluated  and  found  to  be  less  than  1 57c  for  the  following  target  parameters: 

m  =  index  of  refraction  =  1 .5, 
target  volume  <=  1  (pm)3, 

X  =  wavelength  of  light  =  2rt  (pm). 


4  A.  THEORY 

A  brief  outline  of  the  theory  may  be  u-.ef H 1 .  The  outgoing-wave  solution  to  the  Maxwell's  equation  is 


E(r)  =  E0  exp(ikQr)  +^x  Vx  |  d3r'  [(m2  -  1  ),(4^)j  U(r’)  E(r') 
x  exp(ik0'r  -  r‘j)/|r  -  r*|  -(m2  -1)E(r) 


(10) 


:o 


where  L'(  r)  is  one  inside  the  target  and  zero  outside  the  target  to  allow  the  integration  to 
extend  over  all  space.  The  electric  field  is  next  expressed  as 


£.(0  =  Ajj  Eeff>j(r)  (11) 

where  A  is  the  polarization  matrix  for  the  target.  This  matrix  represents  the  effect  of  the  polarization 
charge  induced  within  the  target  by  the  field,  and  it  would  be  expected  to  be  a  function  of  the  internal 
space  of  the  target  and  the  field  variables.  In  the  limit  of  long  wavelengths,  however,  this  matrix  takes  a 
very  simple  form  -  a  constant  diagonal  matrix  for  spheres  and  infinite  cylinders.  Since  we  do  not  know 
w  hat  A  is.  we  assume  a  simple  form  as  the  basis  for  an  iterative  approach.  This  form  is 

A  =  unit  matrix  ,  outside  the  target 

ajE  0  0 

A  =  0  ajE  0  ,  inside  the  target  (12) 

0  0  ayy 

Here,  i  =  3  corresponds  to  the  target  axis  (TE  means  transverse  electric  and  TM  means  transverse 
magnetic).  The  effective  field  is  next  expanded  in  powers  of 

p  =  (m2-  1  )/(4tt)  (13) 

as 

Eeff(r)  =  E0exp(ik0r)  +  In  pn  E(n)eff  (r)  ,  (14) 

and  this  is  substituted  into  equation  ( 10)  using  equations  (11)  and  (12).  Each  term  in  the  series  can 
then  be  determined  from  the  previous  term  starting  with 

E(°)eff  (r)  =  E0  exp(ik0  r ) . 

This  is  the  Shifnn  approach  that  we  have  generalized  for  a  plane  wave  incident  on  a  target  from  an 
arbitrary  direction  and  with  arbitrary  polarization. The  details  of  the  method  along  with  the  first  two 
orders  of  iteration  are  given  in  references  1 2  and  1  3.  For  example,  the  first-order  term  is 


E(1Wr>  =  0<o2  exp(ik0r)/r )  u(k0  r  -  k0)  x[aTE(Eoxt  it  +  Eoyt  jt )+  aTM  E  0zt)]. 


(15) 


where  x  =  rx  r  -  x  ,  and  the  function  u  is  the  Fourier  transform  of  the  pupil  function  for  the 
target: 

u(x)  =  jd'V  U(r)  exp(i  x  r').  (16) 

The  forms  of  u  for  finite  cylinders  and  spheroids  are  given  in  reference  1 3. 

4  B.  THE  AVERAGING  PROCEDURE  FOR  FINITE  CYLINDERS  AND  SPHEROIDS 

We  now  turn  to  the  averaging  procedure  for  cylinders  and  spheroids,  which  is  actually  easier  than 
for  very  long  cylinders  as  discussed  in  section  3.  Refering  to  figure  10  ,  we  recall  the  reference  frame 
(x0  ,  y0  ,  zQ)  ,  where  the  incident  direction  is  z0  and  xQ  is  in  the  scattering  plane.  The  target  frame 

has  zt  along  the  target  axis  and  xt  in  the  plane  containing  zt  and  zQ  (xt .  yt ,  zt )  .The  incident  field 
makes  an  angle  a  with  the  reference  plane  .  and  so  the  components  of  Eq  in  the  target  frame  ,  required 
for  equation  ( 1 5)  are  given  by 


=  E, 


Eoxt 

E0  cos  a 

Eoyt 

—  FI  y[(0) 
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(17) 


The  scattered  w  ave  is  next  obtained  using  the  Shifrin  series  in  the  target  frame  of  reference,  and  it  must 
be  transformed  by  the  inverse  of  the  transformation  used  in  equation  (17)  back  to  the  reference 
frame. These  components  are  denoted  (Esxo,  ESyo>^szo)  Finally,  a  frame  x',  v',  z'  is  define  with  z'  in 
the  scattered  direction  and  x',  z  in  the  same  plane  as  xQ,  zQ  ,  and  the  components  of  the  scattered  field 
are  related  to  this  frame: 


Esx' 

^sxo 

^sy’ 

=  Ryo^o^ 

Esyo 

^sz’ 

^szo 

We  can  now  define  two  intensities: 


(18) 


I  =  (k02  r2  /  Eq2)  [|Esx’|2  +  |Esz.|2P], 

I  =  (k02  r2  /  E02)  |Esy.|2-  (19) 

These  represent  the  intensities  for  the  scattering  of  an  incident  beam  of  light  from  a  target  with  an 
arbitrary  orientation.  For  light  incident  on  a  cloud  of  randomly  oriented  targets,  all  possible  orientauons 
must  be  included.  Thus,  an  average  over  the  orientation  variables  must  be  taken 

lave  =  (1/2rc)  JdyJdo  sin  0  |(y,<t>)  (20) 

4  C.  APPLICATIONS 

Figures  1 1  through  14  show  the  averaged  intensities  lj  and  I  o  versus  0  for  the  scattering  of 
linearly  polarized  light  from  short  cylinders  and  spheroids. The  incident  polarizarion  is  in  the  reference 
plane,  and  the  respective  aspect  ratios  are  10,  2,  1,  and  0.1.  The  targets  have  a  volume  1  (|im)-  .  This 
volume  is  chosen  so  that  the  first  order  in  the  iterative  method  accounts  for  at  least  85%  of  the 
scattering  effect  as  show'n  in  reference  13. 

The  amplitudes  for  the  sphere  are  shown  as  dashed  lines,  while  the  amplitudes  for  the  finite  cylinder 
are  shown  as  solid  lines  with  dots.  The  results  for  spheroids  are  shown  as  solid  lines  with  crosses. 
The  results  for  a  sphere  are  shown  on  each  graph  as  a  reference. 

For  the  average  intensity  I ,  the  results  for  a  spheroid  and  a  finite  cylinder  of  the  same  volume  and 
aspect  ratio  are  quite  similar  except  for  the  minimum  point.  The  long  spheroid  behaves  like  a  cylinder 
while  the  short  spheroid  behaves  like  a  sphere.  On  the  other  hand,  intensities  I  for  the  finite  cylinder 
and  spheroid  differ  markedly  for  the  shorter  aspect  ratios  (the  values  for  the  spheroid  are  zero  for  the 
aspect  ratio  of  one,  and  they  are  too  small  to  appear  on  the  graph  for  the  aspect  ratio  0.1. Therefore, 
spheroids  and  finite  cylinders  of  the  same  material  and  volume  create  similar  scattering  patterns  if  the 
plane  of  linear  polarization  is  not  changed  by  the  scattering,  if  the  orientations  are  random  ,  and  if  an 
average  is  taken.  If  the  plane  of  polarization  is  changed  by  the  scattering,  even  the  averaging  process 
does  not  obliterate  differences  between  the  two  shapes.  Reference  6  details  the  calculations  and 
results. 

5.  SCATTERING  FROM  TARGETS  FORMED  FROM  AN  ASSEMBLY  OF  CIRCULAR  DISKS 

In  this  section  we  will  describe  a  new  model  for  computing  the  scattering  characteristics  of  dielectric 
objects  whose  shapes  may  be  quite  complex.  The  model  assumptions  are:  1 )  the  target  has  a  central 
axis  that  may  be  curved,  and  the  cross  section  at  each  point  on  this  axis  is  circular,  and  2)  the 
wavelength  is  larger  than  the  target  dimension  that  can  produce  standing  waves.  This  second 
assumption  can  be  relaxed  if  a  better  polarization  matrix  is  used  to  initiate  the  calculation  as  will  be 
discussed  later.This  calculational  model  is  then  applied  to  two  cases  whose  scattering  characteristics  are 
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Figure  12.  The  same  as  in  figure  1 1  for  an  aspect  ratio  of  2. 
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v, .  n  t'rom  other  computational  schemes,  tirst  a  finite  dielectric  c>  Under  and  secondly  a  diekxim 

spheroid  (whose  scattering  characteristics  can  be  computed  exactly ).  The  technique  is  then  applied  '■> 
the  case  of  a  dielectric  helix.  In  this  application  the  Circular  Intensity  Differential  Scattering  (C1DS  i  i . 
computed  and  the  results  are  compared  with  the  calculations  of  Bustamente  et.  al.  who  computed  C1DS 
for  a  helix  made  of  a  string  of  dipoles.  The  similarities  between  the  two  calculations  is  striking.  A 
second  senes  of  computations  involve  the  scattenng  of  linearly  polarized  light  by  a  helix.  The  panicle 
parameters  were  chosen  to  match  the  experimental  panicle  used  by  Dr.  R.T.Wang  who  measured  the 
scattering  of  microwaves  by  a  dielectric  helix.  Once  again  the  comparison  between  the  calculation  and 
the  expenmental  results  are  quite  good,  thereby  lending  further  credibility  to  the  calculational  scheme. 
The  final  application  of  this  section  will  be  a  comparison  of  the  scattenng  of  light  from  a  toroid  and 
oblate  spheroid  of  the  same  volume  and  similar  radius  and  height.  We  study  the  signatures  of  these 
shapes  for  different  wavelengths  of  incident  light. 

We  consider  the  scattenng  of  light  from  non-magnetic  dielectric  targets  with  an  index  of  refraction  m. 
Maxwell's  equations  lead  to  an  integral  equation  for  the  final-state  electnc  field  satisfying  the  outgoing 
spherical  wave  boundary  condition  1.6,11,1  — la. 

Ej(r)  =  Einc  j(r)  +  a  1  dV  Djj  G(r/)  Ej(f).  (21) 

where  the  integration  is  over  the  volume  of  the  dielectric  target.  The  elements  of  this  integral  equation 
.ire 

Emc(r)  =  E0  e*P(|ko  r)  =  incident  plane  wave, 

Djj  =  d/dXj  d/dXj  +  k02  S,j  ,  (22) 

G(r,f)  =  exp(ik0|r-r' |  )/|r  -  r*  |  , 

a  =  (m^  -  1  )/(4rt), 

where  the  incident  plane  wave  is  assumed  to  be  harmonic  with  the  time  factor  e.xp(ia)t)  ,  and  where 
k0=2tt/X  , 

In  our  previous  papers,  we  used  the  Shifrin  technique  which  takes  as  the  basis  tor  an  iterative  solution 
the  replacement 

Ek(r')  . >  (Ajk^electrostatic  Ejnc^(r)  ,  (23) 

with  the  assumption  that  the  matrix  A  is  constant  over  the  entire  volume  of  the  target.  This  allows  the 
anal  . tic  evaluation  of  the  first-order  effect  through  the  use  of  Fourier  transforms.  '1  he  technique  is 
shown  to  work  well  for  dielectric  spheres,  finite  cylinders,  and  spheroids  provided  the  diameter  of  the 


target  is  smaller  than  the  wavelength  of  the  incident  light.  For  example,  for  a  long  dielectric  cylinder, 
we  have 

ajE  0  0 

(A)electrostatic  =  0  aTE  0  (24) 

0  0  aiM 

with  aTE  =  2/(m^+l)  and  a-^  =  1  . 

However,  for  more  complex  target  shapes,  the  replacement  (23)  cannot  be  expected  to  hold  as  A 
should  be  expected  to  vary  over  the  target.  We  are  therefore  proposing  a  model  that  allows  for  such  a 
variation  while  retaining  in  part  the  Fourier  transform  technique  If  the  target  has  a  circular  cross  section 
at  every  point  along  its  central  axis,  we  can  divide  the  target  into  disks  which  are  perpendicular  to  the 
axis  and  of  infinitesimal  thickness.  We  then  assume  that  the  polarization  matrix  is  constant  within  each 
disk.  However,  the  matrix  is  permitted  to  vary  from  disk  to  disk.  Thus,  the  integration  in  equation 
(21)  is  divided  into  disks  and 

Ej(r)  =  Einc  j(r)  +  a  I  J  dV*  Dy  G(r.r')  Ajk<  Ek(f)  (25) 

where  the  sum  is  over  the  disks  that  compose  the  target,  and  the  integral  is  over  the  1th  disk  in  which  A! 
is  constant. 

An  expansion  is  then  based  on  the  replacement 

Ajk'  .  Ajk'  Ejnc,k(r)  (26) 

where  A^  is  a  constant  matrix  tailored  to  the  application  as  the  second  model  assumption.  The 
expansion  is  in  powers  of  a  : 

E(r)  =  Einc(r)  +  I  an  E<n)  (r).  (27) 

5  A.  THE  FIRST  ORDER  CONTRIBUTION. 

Putting  the  expansion  (27)  into  equation  (25)  gives  the  first-order  term: 

Ei(1)(r>  -  I  I  dV  [Djj  G(r.r')  ]  Ajk' E,nc  k(r') 

=  Il'(r)  .  (28) 

where  the  summation  is  over  the  disks.  The  disk  integral  I*  can  be  evaluated  by  means  of  Fourier 
transforms  if  it  is  first  converted  to  an  integral  over  all  space  as 


IjV)  =  /  dV*  [Djj  G(r.r*)  ]  Ajk' E0ik  exp(i  ko  0  U(r'), 


(29) 


where  the  pupil  function  U(r')  is  unity  iff  is  within  the  target  and  zero  if  r'  is  outside  the  target. 
Defining  the  Fourier  transforms  by 

G(r.r’)  =  1/(2rt)3  j  d3p  g(p)  exp(i  p  (r-f)  ), 

Djj  G(r,f)  =  1/(2jt)3  J  d3p  g(p)  [k02  Sjj  -  pj  pj]  exp(i  p  (r  -  f)  ),  and 

U(f)=  1/(2tt:)3  J  d3q  u(q)  exp(  iq  r* ),  (30) 

we  evaluate 

g(p)  =  9(P)  =  47t/(p2  -  kQ2),  (31) 

and 

u!(q)  =2?t  a|2  dz|  J  -j  (q  a|)/(q  at)  exp(iqh|).  (32) 

Equation  (10)  is  the  Fourier  transform  of  the  pupil  function  for  a  disk  of  radius  aj  and  thickness  dzp 
The  vector  q  has  the  component  q  perpendicular  to  the  axis  of  the  disk,  and  h  is  the  vector  from  the 
center  of  the  observers  coordinate  system  to  the  center  of  the  1th  disk. 

The  assumpion  that  the  matrix  is  constant  over  the  disk  allows  us  to  integrate  over  the  r’  space 
leaving,  after  application  of  the  convolution  theorem 

lj'(r)  =  1/(2rr)3  J  d3p  u*(p  -k0)g(p)(ko2  Ajk'  -  pj  pj  Ajk>)  Eo  k  exp(i  p  r ).  (33) 

We  now  make  the  far-field  approximation  to  complete  the  integration,  choosing  the  z  direction  of 
p-space  to  be  along  the  unit  vector  r  -  directed  from  the  center  of  the  target  to  the  field  point: 

Ij1  =  1/(27t)3  j  p2dp  j  dq  jdcos(<J> )  exp(iprcos(p  ))  u!(p  -  k0)  g(p) 

A|k  ‘  Pi  Pj  Ajk')  Eo,k  • 

Integration  by  parts,  ignoring  terms  of  order  1/r2  and  higher  gives 

I,1  =  k02  {exp(i  k0  r)/r }  u*(k0  r  - 1^)  [Aikj  -  r,  rj  Ajk*  j  E0  k  ,  (34) 

where  rj,  Tj,  and  rj  are  orthogonal  unit  vectors  with  ^  aligned  along  the  axis  of  the  l1*1  disk. 
Substitution  of  equation  (34)  into  the  expression  for  the  first-order  field,  equation  (28),  gives 


J 


E,(1)(r)  =  2rc  { jds  a2(s)  J-|  (q  (s)  a(s))/(q  (s)  a(s) )  exp(i  q(s)  h(s)  ) 

ko2  (Aik(s)  *  ri  rj  Ajk(s)  )  Eo.k  ) }  e*P('  ko  r)/r-  (35) 

The  vector  q  is  the  transfer  of  momentum  during  scattering,  kGtr  -  kG)  .  with  the  unit  vector  kG 
pointing  along  the  direction  of  incidence.  In  equation  (35)  the  sum  over  1  has  been  replaced  by  an 
integration  along  the  main  axis  of  the  target,  denoted  by  s.  We  will  see  in  the  applications,  that  the 
frame  of  reference  to  be  used  is  suited  to  the  1th  disk  with  the  disk-z  axis  along  the  axis  ot  the  disk.  In 
tins  frame,  A  will  be  assumed  constant  and  diagonal. 

5  B  THE  SECOND  ORDER  CONTRIBUTION. 

Returning  to  the  expansion  of  the  final-state  electric  field,  equation  (25). the  second -order  term  is 


E,(2)(r)  =  I  JdV'  [DjjG(r,r')]  Ajk*  Ek^)  (r')  U|(r'). 

E(  ‘  ’  is  the  first-order  term  within  the  target: 

Ek0)(r)  =  I  lk'(0 

»  l/(2n:)3  I  J  d3pu'(p  -  kg)  g(p) 


(ko  Akm  '  Pk  Pn  Anm  )  Eo,m  exP('Pr) 


introduction  of  the  Fourier  transforms  in  eqs.  (30)  allows  the  integration  of  r'-space  leaving 

e/2)  (r)  =  1/(2it)3  I  I  Jd3p’  Jd3p  g(p')  [k02  8a  -  p’,  p\  ] 


xu  (p'  -  p)  u  (p  -  kQ)  g(p) 
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We  next  make  the  far- field  approximation  which  allows  the  lnteciam-n  •  we-  p  ace  using  integration 
1  >v  parts  iitnorme  terms  of  the  order  of  I/m.  the  result  is 


E32)  (r)  =  exp(i  k0  r)/(4^r)  I  I  Jd3p  g'D>  ;k0“:  o. 

xu 1  ( p  -  k0  r  )  u!(k0  r-k0)  A]k' 
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Figure  15.  The  reference  frame  (0)  and  the  target  i runic  (ti.  Light  is  incident  along  the  z()  axis  and 
scattered  along  r.  The  symmetry  axis  of  the  target  is  along  /  and  this  axis  is  oriented  by  the  polar  angle 

0t  and  the  azimuthal  angle  yx. 
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The  summations  will  become  integrations  over  the  symmetry  axis  of  the  target,  and  so  there  is  a 
five-fold  integration  to  be  performed 

6  .  APPLICATIONS  TO  TARGETS  WITH  A  LINEAR  SYMMETRY  AXIS 
6  A,  THE  FINITE  DIELECTRIC  CYLINDER 

The  results  of  this  application  without  the  decomposition  into  small  disks  have  already  been 
presented  in  Refs.  6  and  13  .  Although  there  is  no  need  to  divide  this  target  into  infinitesimal  disks 
when  the  polarization  matrix  is  assumed  to  be  constant  within  the  cylinder  (Shifrin  approximation),  it  is 
nevertheless  instructive  to  use  the  cylinder  to  show  how  the  disk  technique  works.  On  the  other  hand, 
this  decomposition  would  be  necessary  if  we  wanted  to  apply  the  results  of  reference  13  which 
demonstrated  that  A  does  indeed  vary  over  the  cylinder  volume  for  short  cylinders. 

The  cylinder  geometry  is  shown  in  Fig.  15.  The  reference  frame  (x0,y0,z0)  is  defined  with  zQ 
along  the  direction  of  the  incident  light,  and  the  Xq  axis  is  chosen  such  that  the  scattered  direction  r.  lies 
in  the  X()-Z()  plane.  The  target  frame  txt,yt.zt)  is  defined  with  zt  along  the  symmetry  axis  of  the 
cylinder,  and  the  Xt  axis  is  chosen  so  that  the  direction  of  incidence,  zq  lies  in  the  x,-z.  plane.  If  the 
polar  and  azimuthal  angles  of  zt  are  denoted  0t  and  yt  ,  respectively,  the  reference  and  target  frames 
are  related  as 

Xt,i  =  T0->t,ij  *o.j  (40) 

with  the  transformation  matrix  given  by 

cos  <t>j  cos  cos  Of  sm  y^  -  sin 
T0->  t  =  'sin  ‘ft  cos  yt  0 

sm  c>t  cos  yj  sm  Oj  sm  yt  cos  0{  (41 ) 

In  the  target  frame  of  reference  the  polarization  matrix  is  both  constant  and  diagonal: 

a-|-£  0  0 

A  =  0  ayj=  0  (42) 

0  0  aTM 

For  a  long  cylinder  (aspect  ratio  =  length/diameter  >  20),  the  matrix  elements  in  the  long  wavelength 
limit  are  given  in  equation  (4)  For  shorter  cylinders,  reference  13  gives  values  obtained  in  the 
long-wavelength  limit  where  the  values  tor  disks  nearer  the  edges  differ  significantly  from  those  nearer 
the  center 


The  first  order  contribution  to  the  final-state  electric  held  has  its  components  first  expressed  m  ti  e 
target  frame  to  take  advantage  of  the  form  in  equation  (22 ): 

E^(r)  =  2rt  a2  k02  J-|(q  a)/(q  a)  exp(i  k0r)/r 

x  [ajE^oxt 't  +  EoytJt)  +  aTM  Eozt  kt  transverse 

xj  dz  exp(i  q  z  )  (43) 

Here,  "transverse”  denotes  that  portion  of  the  vector  transverse  to  the  scattered  direction  of  light: 


^transverse  -  X  -  X  r  r 


Moreover, 


q  =  |k0  (r-k0)x  kt  |, 

q  =  kQ  (r  -  ko)  i  (44) 

where  kt  is  a  unit  vector  along  the  symmetry  axis  of  the  cylinder,  a  is  the  radius  of  the  cylinder,  and  h 
is  the  height  of  the  cylinder. 

The  incident  wave  is  assumed  to  be  linearly  polarized  in  this  application: 

Eo  =  Eoxo  'o  +  Eozo  Jo  (45) 

in  the  reference  frame,  and  EQ  makes  an  angle  a  with  the  xD  axis.  The  components  of  EQ  used  in 
eqation  43  are  in  the  target  frame,  thus  the  transformation  (20)  is  used.  Evaluation  of  the  integral  in 
equation  43  gives  the  result  derived  in  reference  13: 

Ej(^(r)  =  2:ta2  h  J i  (q  a)/(q  a)  exp(i  kQ  r)/r  kQ2 
IaTE<Eoxt 't  +  Eoyo  it)  +  aTM  Eozt  ktl 

sin(q  h/2)/(q  h/2).  (46) 

As  the  height  of  the  cylinder  increases,  the  last  term  will  give  a  significant  contribution  only  when  q  = 
0.  We  show  in  reference  6  that  this  gives  rise  to  the  scattering  being  restricted  to  a  volume  near  the 
surface  of  a  cone  whose  symmetry  axis  lies  along  the  cylinder  axis  and  whose  apex  halt-angle  is  o,, 
The  scattering  angle,  between  knand  r .  is  denoted  by  20o. 

This  derivation  shows  that  the  division  of  the  target  into  disks  and  the  subsequent  coherent  addition  of 


these  effects  leads  to  the  same  first-order  result  as  that  derived  for  the  scattering  of  light  from  a 
cylindrical  target.  We  should  also  emphasize  that  this  new  approach  is  necessary  if  the  polarization 
matrix  vanes  along  the  cylinder  axis  with  cylindrical  symmetry  as  was  shown  to  be  the  case  in  the  same 
reference.  A  calculation  is  now  performed  for  a  cylinder  of  volume  4.8  (|itn)3  and  aspect  ratio  5  for 
incident  light  of  wavelength  1  =  2n  microns.  The  incident  light  is  linearly  polarized  in  the  x0-z0 
plane.  These  conditions  insure  that  the  first-order  calculation  is  dominant  and  accurate  as  shown  in 
reference  13. The  polarization  matrix  elements  vary  over  the  cylinder  along  the  symmetry  axis  x  as 

aTNl  =  *  '  0-032  (x  -  ar/2)“. 
aTE  =  0.722  +  0.016  (x  -  ar/2)2  , 

where  ar=5.  These  equations  match  the  numerical  results  contained  in  figure  7  of  reference  13  The 
intensity  hi  is  calculated,  where  both  the  incident  and  scattered  light  are  polarized  in  the  xQ-z0  plane, 
and  the  results  are  show-n  here  in  figure  16. 

o  B.  THE  DIELECTRIC  SPHEROID 

We  next  turn  to  an  application  where  an  exact  calculation  is  available  for  comparison  14  We  have  also 
performed  this  calculation  using  the  Shifrin  method  since  the  electrostatic  polarization  is  known.  *3  The 
disk  method  is  therefore  applied  to  the  dielectric  spheroid  to  establish  its  validity,  and  the  range  of 
v  alidity  of  the  first-order  calculation. 

The  geometry  is  the  same  as  used  for  the  finite  dielectric  cylinder.  The  first-order  calculation  is 
similar  to  that  for  the  cylinder  except  that  the  radius  of  the  disks  vary  along  the  symmetry  axis.  The 
first-order  contribution  to  the  final-state  electric  field  is 

=  27tk0^  exp(i  kQ  r)/r 

laTE(Eoxt't  +Eoytlt)  +  aTM  E  ozo  ^transverse 
j  dza2(z)J1(q  a(z))/(q  a(z))  exp(i  q  z).  (47) 

a  here 

a(z)  =  a/b  [b2-z2] 1 /2  ,  (48) 

tor  both  oblate  and  prolate  spheroids.  Here,  a  is  the  radius  of  the  centra!  circular  cross  section,  and  b 
the  half-height  of  the  spheroid.  The  polarization  matrix  elements  aTp  and  a  j  are  w  ell  know  n  in  the 
ring  wavelength  limit  and  quoted  in  reference  1 3  for  prolate  and  oblate  spheroids. 

The  intensites  1 1  ]  and  In  are  defined  as 

'll  =  (kor/Eo  )2  iEsc|2ff  =  90  • 


where  a'  is  the  angle  of  incident  linear  polarization  measured  with  respect  to  the  x0-z0  plane.  The 
first-order  electnc  field  is  calculated  using  equation  (27)  and  compared  with  the  exact  results  of  Asano 
and  Yamamoto,  reference  14,  using  a  spheroid  w  ith  dielectric  constant  m  =  1.33.  volume  V  =  4.X1 
pma  and  the  aspect  ratio  ( height/diameter)  5.  The  wavelength  of  the  incident  light  is  >.  =  2k 
microns.  The  result  for  In  is  shown  in  figure  16  along  with  that  for  the  comparable  cylinder 
discussed  above.  It  is  first  noted  that  this  result  for  the  spheroid  is  in  excellent  agreement  with 
reference  14-  as  it  was  in  reference  13.  In  reference  i 3.  the  Shifnn  method  was  used  to  calculate  the 
first  and  second  order  contributions  to  the  final-state  electric  field.  It  was  found  that  the  second  -order 
contribution  is  not  important  for  this  target  size,  wavelength,  and  index  of  refraction. 

More  important  for  the  purpose  of  this  report,  this  calculation  show  s  that  the  decomposition  of  the 
target  into  disks  is  accurate  when  the  disk  radii  vary.  This  allows  one  to  approach  applications  that 
cannot  be  done  by  using  the  Shifrin  method,  nor  any  approach  other  than  a  brute-force  evaluation  of  the 
starting  integral  equation,  equation  (25).  But  even  for  these,  one  needs  to  know  the  polarization  matrix. 

One  also  notes  in  figure  17  that  the  cylinder  and  the  spheroid  have  distinct  scattering  signatures  for 
this  choice  of  scattering  parameters. 

7.  APPLICATION  TO  DIELECTRIC  TARGETS  WITH  A  WANDERING  SYMMETRY  AXIS 
7  A.  THE  DIELECTRIC  HELIX 

We  now  treat  an  extended  target  that  has  locally  circular  cross  sections  but  whose  axis  is  a  helix.  The 
geometry  is  complicated  by  this  extension,  and  three  frames  of  reference  are  needed  to  describe  the 
scattering  process.  The  reference  frame  (xQ,  yQ,  zQ)  is  as  described  above.  The  target  frame  (xt,  yt,  zt) 
takes  zt  along  the  helical  axis.  A  third  frame  of  reference  (Xq  ,  y^  ,  z0  )  called  the  disk  frame  is  placed 
on  a  disk  of  infinitesimal  thickness  with  z^  through  the  center  of  the  disk  and  perpendicular  to  it.  This 
geometry  is  shown  in  the  three  pans  of  figure  17  . 

The  axis  of  the  helix  is  traced  by  the  vector  h  defined  in  the  target  frame  as 

h(o)  =  R  cos  o  ^  +  R  sin  o  it  +  P/(2tc)  o  kt  (50) 

Here.  R  is  the  radius  of  the  helix,  P  is  the  pitch,  and  0  is  the  angular  variable  that  is  shown  in  the 
figures  18.  The  +  and  -  signs  define  right-  and  left-handed  helicity.  respectively.  Note  that  0  runs  from 
0  at  one  end  of  the  helix  (  at  the  origin)  to  2itn  (  n=  number  of  turns)  at  the  other  end. 

The  vector  h ( 0 )  marks  the  position  of  the  center  of  a  disk.  The  disk  z  axis  is  along  the  differential 
vector 

dz  =  db/do  do 

=  {-R(sin  o  if  -  cos  0  j| )  +  P/(2n:)  kt  }  do.  (51) 


o 


Figure  17.  A)  The  reference  and  target  frames  of  reference  shown  for  a  helical  target.  B)  A  detailed 
view  of  the  helix  in  the  target  frame  of  reference.  An  infinitesimal  disk  is  shown  located  by  the  vector 
h(0)  with  azimuthal  angle  Q.  C)  A  detailed  view  of  a  disk  in  the  disk  frame  of  reference.  The  disk  radius 
is  denoted  a. 
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The  magnitude  of  dz^  is  independent  of  d  and  given  by 

dz  =  [R2  +  (P/27t)2]1/2d<{>.  (52) 

The  unit  vectors  for  the  disk  frame  are  chosen  as 

k  =  {-R  (sin  d  it  -  cos  <|>  jf )  +  (P/2tc)  kt }  /S  , 

i  =  cos  d  if  +  sin  d  jf  , 

j  =  k  x  i  =  {+  (P/2rt)  (cos  d  it  •  sin  d  if)  -  R  kt  }  /S, 

S  =  [R2  +  (P/2tc)2  ]1/2.  (53) 

The  transformation  from  the  reference  to  the  disk  frame  is  achieved  through  the  orthogonal 
transformation 

xd,i  =  ("*"t->  d  )ij  C^o ->  t  )jk  xo,k  =  (To  ->d  )ik  xo,k-  (54) 

The  transformation  matrix  is  given  in  the  appendix  of  reference  Note  that  this  transformation  is  unitary, 

Td  ->  0  =  T0  ->  d 

Our  model  assumes  that  the  polarization  matrix  A,  appearing  in  equation  (1),  is  constant  and  diagonal 
in  the  disk  frame  of  reference,  with  the  x^  and  y^  matrix  elements  equal  due  to  the  cylindrical 
symmetry.  To  take  advantage  of  this  assumption,  each  term  in  the  integrand  of  equation  (15)  must  be 
expressed  in  the  target  frame  and  then  transformed  to  the  reference  frame.  The  first-order  contribution 
to  the  field  is  therefore  given  by: 

Ej(1)(r)  =  2xexp(ik0r)/r  k02  /  S  dd  a2  Jf  (q  (d)  a)/(q(d)  a) 

exp(iq(d)  h(d)  )  (T  <j>  — >  o ) ij  rj  ^k  ^kl)  ^o,l(d)  Idisk’  (55) 

where  n  is  the  number  of  turns  in  the  helix.  The  transformation  matrix  Tf  _>  0  is  introduced  in 
equation  54  and  specified  in  the  appendix  of  reference  Here,  a  is  the  radius  of  each  disk  (wire  radius) 
and  q  =  kQ  (r  -  kQ),  r  is  a  unit  vector  directed  from  the  origin  of  the  target  frame  to  the  field  point  (the 
far-field  approximation  is  used)  and  kQ  is  a  unit  vector  in  the  direction  of  the  incident  light.  Thus,  it  is 
necessary  to  continuously  transform  from  the  disk  to  the  reference  frame  as  the  integration  over  the 
heiix  is  performed.  Therefore,  the  scattering  effects  from  the  disks  that  compose  the  target  are  added 


coherently,  their  phases  being  matched  by  the  exponential  factor  in  the  integrand  of  (55).  Moreover, 
the  vector  q  appearing  in  the  Bessel  function  is  the  magnitude  of  the  vector  q  perpendicular  to  the  local 
disk  axis.  We  get  this  quantity  by  using  the  transformation  of  equation  (54): 


qx 

ko  sin(2  <>o) 

qy 

=  T 

0 

k0  cos(2  q0) 

where  2  0o  is  the  scattering  angle.  It  follows  that 


(56) 


q  (4>)=  [(q><t,  )2  +  (qjq  )2  ]1/2 
=  {(cos  qq^  +sin0qyt  )2  + 

(-R  qzt  +  (P/2n)(cos  q  q^  -  sin  q  qxt  )2  /S2  }1/2  (57) 


The  phase  term  in  equation  (35)  has  the  exponent 


q  h  =  R  cos  q  qx  +  R  sin  q  qy  +  (P/2n)  q  qz 


(58) 


The  last  factor  in  equation  (35)  contains  the  components  of  the  incident  field  in  the  disk  frame; 
therefore,  the  incident  electric  field,  normal  to  k0  ,  must  be  transformed  as 
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(59) 


7  B.  CALCULATION  OF  THE  CIDS  FOR  THE  DIELECTRIC  HELIX 
The  scattering  of  light  from  very  thin  helical  targets  has  been  discussed  by  various  authors  J5-18 
where  the  helix  is  treated  essentially  as  a  line  of  dipoles  set  on  the  helical  axis.  Singham.  Patterson, 
and  Salzman,  reference  19,  have  also  considered  light  scattering  from  chiral  objects  ,  including  the 
helix,  using  a  model  that  descibes  the  three-dimensional  body  as  a  collection  of  spherical  or  ellipsoidal 
dipoles.  The  quantity  of  special  interest  in  these  papers  is  called  the  circular  intensity  differential 
scattering  parameter  (CIDS )  defined  as 


CIDS  =  (l[_-  Ip)/  (lj_+  Ir)  . 


(60) 


Here  Ip_  is  the  scattering  intensity  for  left-circularly  polarized  light 


E0  ~  Eq  ('o  ±  '  Jo)’ 


with  the  upper  and  lower  signs  corresponding  to  right  and  left  circularly  polarized  light,  respectively 
(with  the  time  factor  icot  ).  In  the  series  of  papers  by  Bustamante  and  others,  refs.  16-18  ,  the  CIDS 
parameter  is  evaluated  for  very  thin  wires  for  which,  in  our  notation,  aTE  =  0  and  aTN1  =  1.  In 
reference  16.  a  series  of  calculations  is  performed  for  a  helix  with  n=  20.  R=  1pm  ,  p  =  1  jam,  and  for 
wavelengths  from  1  |im  to  4pm.  The  spiral  axis  is  located  at 

Of  =  90  and  yt  =  90 

We  have  calculated  the  CIDS  ratio  for  the  helix  described  above  for  the  wavelengths  1,2,  4,  and  6 
urn.  It  should  be  emphasized  that  our  model  is  different  than  that  of  Bustamante.  Our  helices  have 
thickness  and  are  dielectrics  of  refractive  index  1.5.  The  polarization  matrix  is  diagonal  and  of  the  form 
of  equation  (22)  with  a-^  =  2/(m~  +  1)  and  aTN!  =  1  in  each  disk  frame  of  reference.  In  general,  we 
find  that  our  results  give  about  the  same  signatures  as  those  of  Bustamante  et  al.,  that  is  the  number, 
size,  and  placement  of  the  lobes  are  about  the  same  as  found  by  them.  However,  we  do  not  get  their 
large  back*  ottering  lobes  .  FigurelB  is  a  polar  plot  of  the  CIDS  ratio,  where  1  =  1  (am.  Figure  19  is 
similar  b  .ah  1  =  2pm.  In  figure  20,  we  attempt  to  match  the  Bustamante  calculation  by  taking 
m=50  v.  rich  makes  aTC  small  (  note  that  the  radius  is  correspondingly  small  so  that  the  condition 
2  tram  A  <  2  is  satisfied  ).  After  doing  this  we  get  a  large  backscanenng  Jobe,  figure  21  is  a  CIDS  plot 
for  1  =  4  pm.  m=1.5  ,  and  only  one  prominant  lobe  is  present,  and  again  this  is  consistent  with 
reference  16.  Figure  22  is  calculated  as  in  figure.  6  but  with  1  =  4  pm  ,  and  a  backscattenng  lobe 
again  appears.  Figure  23  is  a  CIDS  plot  for  a  wavelength  1  =  6  pm.  Here  there  is  one  very  wide  lobe 
with  a  large  backscattenng  contribution  even  though  the  calculation  is  for  a  thin  dielectric. 

It  should  be  pointed  out  that  o  -  results  for  a  thin  dielectric  differ  from  the  results  of  reference  16  even 
when  the  radius  of  the  wire  is  •  a  to  be  very  thin,  a  =  0.000 1  pm.  and,  in  fact.  CIDS  does  not  change 

significantly  for  radii  from  thi  *  e  to  about  0.1  pm.  Thus,  a  thin  dielectric  wire  is  not  the  same  as  a 
tram  of  dipoles  placed  along  t:  al  axis  .  Any  fmue  thickness  seems  to  make  a  difference,  especially 

in  the  back  scattering  range.  '1  suit  seems  to  be  borne  out  by  the  work  of  reference  18  ,  where  it  is 

found  that  tnaxial  poiarizabih  (which  are  related  to  differing  diagonal  matrix  elements  of  the 

polarization  matrix)  significant!,  affect  the  scattering  result,  and  therefore  have  to  be  included  to 


qualitatively  describe  the  scattering  matrix.  The  model  of  reference  lx  .  though  not  treating  the  target  as 
being  a  continuous  body  as  does  our  model,  is  nevertheless  three  dimensional  and  thus  closer  to  our 
treatment  than  the  one  dimensional  model  of  reference  16. 

One  final  calculation  is  made,  based  on  .in  electron  micrograph  of  Spirillum  membranes  presented  in 
’'cteience  lb  These  researchers  use  the  following  scattering  parameters  *<>r  this  application:  /.  -  260 
tiro,  pitch  =  2.6  nm,  radius  of  helix  =  208  run.  and  b  turns.  Our  calculation  uses  a  real  index  of 
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Figure  23.  The  CIDS  parameter  for  the  same  helix  as  in  figure  IS.  The  wavelength  of  the  incident  light 
is  6  pirn,  and  the  polarization  matrix  elements  for  each  disk  are  as  described  in  figure  18. 
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Figure  24.  The  CIDS  parameter  for  a  biologically-  modelled  helix.  The  index  of  refraction  is  taken  to  be 
m  =  1.33.  This  helix  has  6  turns,  and  its  size  parameter  and  wavelength  of  the  incident  light  are  given  in 
the  text.  The  polarization  matrix  for  each  disk  is  as  described  in  figure  18. 


refraction  of  1.33.  and  our  results  are  shown  in  figure  24.  The  results  of  refer  etu.  e  >  '  •• 
application  show  two  lobes  with  almost  the  same  maximum  for  the  ('IDS  parameter.  a;x  ,t  A 
Our  results  indicate  a  thinner  lobe  w  ith  the  same  maximum  value  oecuring  at  about  the  same  scatter 
angle  of  80. 

7  C  .  DIELECTRIC  TOROID. 

The  helix  reduces  to  a  toroid  when  the  pitch  p  =  0  Therefore,  the  evaluation  of  the  electric  :.c.  : 

*> 

proceeds  as  for  the  helix.  For  a  dielectric  toroid  with  volume  V  =  4.8  pnr'  .  the  intensity  ;• 
calculated  for  two  wavelengths.  K  =  n  (am  and  2:r  jam.  This  target  si/.e  insures  that  the  first  order 
contribution  is  accurate  by  virtue  of  our  w  ork  in  reference  13. 

The  scattering  intensity  I->->  is  calculated  and  compared  to  tut  oblate  spheroid  of  the  same  volume  at:.: 
aspect  ratio  1/5.  The  index  of  refraction  of  both  targets  is  m  -  1.33.  and  the  toroid  has  a  =  0.358  L.m. 

R  =  1  9  jj.ni.  The  oblate  spheroid  has  a  radius  I  CO  uni  and  a  height  0.716  jam.  Again,  each  shape  has 
the  same  volume.  The  results  are  shown  in  figures  25  and  26.  Figure  2^  shows  quite  different 
signatures  for  the  two  shapes  when  the  wavelength  is  n  jam.  the  toroid  alone  has  a  minimum  at  about  a 
scattering  angle  of  40  degrees,  and  the  backscattenng  intensity  for  the  toroid  is  appreciably  larger  than 
for  the  spheroid  (note  that  this  is  a  log  plot,  so  the  difference  in  backscattenng  intensities  is  almost  two 
orders  of  magnitude >.  In  figure  26.  with  /.  -  2 rc  pm.  these  differences  are  nearly  eliminated  •  the 
or.lv  real  difference  between  the  patterns  being  the  central  minima  oecuring  at  slightly  different 
scattering  angles. 

8.  SCATTERING  OF  LIGHT  FROM  A  DIELECTRIC  HELIX  -  EXPERIMENT 

Experimental  studies  of  scattenng  by  a  single  particle  are  very  difficult  to  perfonn  in  the  optical 

region,  especially  when  the  particle  size  is  comparable  to  the  wavelength  of  the  incident  radiation. 

The  considerable  advantages  of  using  the  microwave  analog  technique  in  simulating  single  particle 

x  in 

scattering  have  been  described  elsewhere.  We  select  here  some  recent  angular  distribution 
measurement  results  for  a  dielectric  helix  employing  this  technique. 

A  right  handed.  7-turn  helix  of  3.66  cm  outer  diameter  and  of  4.35  cm  axial  length  was  formed  from 
a  0.48  cm  diameter  plexiglass  cylinder,  a  commercially  available  acry  lic  rod.  The  complex  dielectric 
constant  of  this  plastic  material  was  separately  measured  by  a  classical  standing  wave  technique  using  a 
waveguide-slotted-line,  from  which  the  complex  refractive  index  was  found  to  be  m  =  1.626-  0.012  i 
at  the  operating  microwave  wavelength  /„=  3.18  cm. 

The  microwave-unique  compensation  (or  null >  technique  was  employed  in  birth  extinction  (20o  =  0  ) 
and  in  angular  scattenng  (10  <20o  <  170  (intensity  measurements,  and  from  th*  former,  one  can 
also  deduce  the  scattering  intensity  in  the  beam  direction  (20o  =  O'  i.  Presently,  the  null  technique  is  the 
onlv  reliable  methixl  with  which  one  can  discriminate  the  true  scattered  wave  against  the  superimposed, 
vet  unwanted  coherent  background  radiation,  whom  magnitude  may  he  comparable  to  or  even  much 
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future  25.  The  intensity  I-,,  tor  a  dielectric  {•  >io;.i  aiui  an  oblate  spheroid.  I  he  volume  of  each  target  is 
4  H  pm^  and  the  index  of  each  is  m  -  !  I  he  other  si/e  parameters  are  given  in  the  text.  The 
wavelength  of  the  incident  light  is  turn 
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larger  than  the  former. 

A  brief  explanation  of  measurement  procedures  using  this  technique  is  needed.  In  the  absence  of  the 
target  in  the  beam,  one  cancels  the  unwanted  background  at  the  output  port  of  a  mixer  bv  properlv 
adjusting  both  the  amplitude  and  the  phase  of  a  reference  wave  which  is  piped  through  a  separate 
waveguide  from  the  microwave  source  to  an  input  port  of  the  mixer,  while  the  background  wave  via 
free  space  is  being  fed  to  the  other  mixer  input  port.  As  the  target  is  hoisted  into  the  beam,  the 
off-balance  appearing  in  the  mixer  port  is  directly  proportional  to  the  true  scattered  wave.  The  target  is 
thereby  allowed  to  go  through  a  preselected  sequence  of  orientations  in  the  beam,  and  the 
onentation-dependent  data  are  recorded.  Immediately  after  the  run,  the  target  is  replaced  by  a  standard 
target  of  known  scattering  magnitudes,  and  the  measured  scattering  of  this  standard  provides  the 
absolute  magnitude  calibration  of  the  entire  run. Therefore,  it  follows  that  the  accuracy  of  the  scattering 
experiment  hinges  most  critically  on  the  stability  of  cancellation,  i.e..  on  how  long  the  initiallv 
established  null  condition  is  maintained.  Any  disturbance  affecting  the  initial  cancellation  will  apperar 
as  drift  or  noise  in  the  detected  signal,  which  is  virtually  impossible  to  correct  after  the  measurement  is 
made.  It  also  follows  that  the  measurement  accuracy  degrades  as  the  scattered  signal  level  becomes 
lower  in  comparison  to  the  drifting  background  or  reference  wave. 

The  measured  scattering  intensities  (in  absolute  magnitude)  are  denoted  I]  j  and  l->s  .  The  intensitv 
I  [  j  is  the  case  when  both  the  receiving  and  transmitting  antennas  are  polarized  vertically.  I->->  is  for 
both  antennnas  polarized  horizontally  The  horizontal  scattering  plane  corresponds  to  the  x0-z0  plane  in 
figure  IT  Inspection  of  the  experimental  profiles  1 1  j  and  bs  (figures  28.29,30)  versus  the  scattering 
angle  2o()  immediately  shows  that  the  scattering  by  the  helical  structure  is  hightly  sensitive  to  the 
orientation  of  the  helix  in  the  beam;  and  also  that  the  intensity  may  vary  as  much  as  four  orders  of 
magnitude  (  figure  28)  as  the  scattering  angle  changes.  Because  of  the  rather  small  size  of  the  helix 
compared  to  the  large  antenna-target  separations  (  10.8  m  between  the  transmitter  and  target;  5.2  m 
between  receiver  and  target  ),  and  hence  the  low  scattered  signal  levels,  significant  measurement  errors 
mas  exist  due  to  the  afore-mentioned  compensation  stability.  Roughly  estimated,  the  cumulative 
maximum  errors  would  be  about  5T  if  the  absolute  magnitudes  of  I  j  j  and  I  ->  ->  are  larger  than  10,  up 

approximately  20  D  for  those  between  10*'  and  10* .  and  may  amount  to  as  much  as  100  O  for  those 
magnitudes  having  less  than  10'*.  Also,  the  system  gam  changes  during  the  lengthy  measurement  runs 
mas  have  contributed  additional  absolute-magnitude-calibration  errors.  Roughly  estimated,  the 
cumulative  maximum  errors  would  be  about  5T  if  the  absolute  magnitudes  of  I  j  j  and  In  are  larger 
than  10*.  up  to  approximately  20  Of  for  those  between  10(l  and  10*.  and  may  amount  to  as  much  as 
10<r;  tor  those  absolute  magnitudes  having  less  than  10'*. 

■»  COMPARISON  OF  THEORY  AND  EXPERIMENT 

I  he  senes  of  measurements  on  the  scattering  of  microwaves  from  a  dielect  ; :  spiral  target  have  been 
pertormed.  a>  described  alios e.  and  theoretical  calculations  are  compared  ssith  the  intensities  I  j  j  and 


Figure  27.  a>  The  reference  frame,  u here  i>  m  t!ie  direction  of  incident  light,  r  is  the  direction  ot 
scattered  light,  and  x„  is  :n  the  scattering  plane  !'• -rated  ‘w  /  r  I  he  axis  ol  the  helix  is  along  zt 
whose  orientation  ancles  are  o,  .aid  1  ne  : -  cate:.:  e.emru  t.ei.l  >n  is  L„  ativi  it  makes  an  angle 


(iwith  the  scatter. me  pi. me  The 


hi  T^e  helix  show  ■ 

■ ,  Detailed  d:  - 


hi o  .  1  .  .a 


s’ 

'T 

e 


Logf^ « (90,0)] 


Scattering  Angle 


Scattering  Angle 


Figure  30.  Experimental  scattering  profiles  I ,  j (90,90),  1:2(9(),90)  versus  scattering  angle  2  Q( 
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i--  tor  w  inch  the  plane  of  polarization  is  no;  changed  bv  the  scattering 
1'he  range  of  parameters  tor  the  helix,  user!  in  these  experiment-  ate  :::  t!:e  ranee  •>;  aivdc 
integral  method  for  calculating  scattering  intensities.  The  intcu-iiic-  cumulated  are  lj  •. 
directions  of  initial  and  final  linear  polarization  are  perpendicular  to  the  scattering  plane.  and  1  - 
the  directions  of  linear  polarization  are  in  the  scattering  plane  i  the  scattering  plane  i-  fottix 
directions  of  incidence  and  scattering  -  the  xt>- z0  plane  of  the  theoretical  disctissa  <r.  • 

The  experimental  results  for  IjjtO.O).  IssiO.Oi  are  show  n  in  figures  2V  tire  result-  tor  1 
and  IssiliO.()i  are  shown  in  figures  29 .  and  the  results  for  Iioui).o;ii  and  I-.mOO.9Ui  arc  - 
f. cures  30.  }  lere.  for  example  Ij  jt90.0i  corresponds  to  the  heitx  orientation  o.  on  and 
o  .  Note  that  the  values  never  drop  to  below  10"*-.  even  though  the  theor\  predicts  mamma  :h 
smaller.  We  account  for  this  by  noting  that  the  real  target  is  not  perfect!)  shaped,  and  hence 
sharp  minima  produced  by  an  ideal  spiral  would  not  be  expected.  Moreover,  the  r.otse  level  m. 
the  experiment  prohibits  meaningful  measurements  three  orders  of  magnitude  below  maximal 
The  positions  of  the  minima  are  of  greatest  importance  as  these  are  produced  bv  the  spiral  c 
and  we  summarized  these  minima  for  the  experiments  from  figures  2X  as  best  we  can  as 


for  futures.  29  as 


and  for  futures.  30  as 


l22(0.0)  :  50n  75',  1 00 7 .  1 40 7 
1^10.0)  :  40 7 .  707  .  1  35 7 . 

I22(90.0):  50°,  70\  90  ,  1  50  ; 
I-,  1(90.0):  50'.  70 : .  1057  140 

li  i  (90,90):  457,  115  ; 
l22(90,90):  557.  907  125m 


f.g. ires  '!  through  3o  contain  the  intensities  ij  ](0.0i  InjOjli.  1]  |  Pin. to.  |-,m90.0'  .  I- 
,md.  I ^ si 90.90 1  and  the  comparison  of  the  theory  and  experiment.  The  experimental  points 
a>  solid  diamonds  while  the  theoretical  results  are  plotted  as  a  line. 

The  comparison  between  theory  and  experiment  is  evident  from  the  figures  m  uhnh  ti  e  • 
results  are  normalized  with  the  experiment  at  the  smallest  scattering  angle  t  10  i.  We  can  s.j 
this  comparison  as  follow s: 

1 1 1  Reasonably  good  agreement  between  theory  and  experiment  tor  forward,  and  back  -cat:, 
forward  to  back-scattering  ratios  tire  particularly  well  established  bv  the  theory  tor  the  : 
IssiO.Oi.  I,  ,190.0).  and  1 22* <v>( )*t>< ’ 1 

'  2  i  ( iooti  agreement  between  theoretical  and  experimental  values  tor  the  I,, 

•  taut. .re  of  the  intensity  patterns  for  1|  j  (0.0)  and  I  j  ]<  90,90/  is  p.irttcularlv  stnkmc.  and 
adeemed!  > veurs  where  the  both  theoretical  and  experimental  vah.e  are  within  two  outers  ■ 

•  g  tiie  forward  -cat  ten  tig  values.  In  general,  agreeemem  i-  hot  w !  ere  t; -•  :!  eorem.  .d  •.  ,f  . 
de--end  below  about  two  or  three  orders  of  magnitude  ot  the  forward,  values.  Pus  nxil:  ..a 


■'  /■. >'c.' s 


since  the  experimental  values  have  large  uncertainties  for  such  small  intensities  as  described  above. 

(3)  The  positions  of  most  of  the  experimental  minima  generally  agree  well  with  the  theory. 

One  concludes  that  the  theory  presented  here  explains  the  general  features  of  the  spiral  experiment 
Since  the  theory  can  be  applied  to  a  wide  range  of  particle  shapes,  this  technique  provides  a  method  b\ 
which  such  target  shapes  can  be  studied  without  resorting  to  time-  consuming  numerical  procedures. 
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Figure  33.  A  comparison  betvsen  theory  and  experiment  tor  the  log  t  I, 
angle  2  Oq  for  the  orientation  Qt  =  90  and  y(=  0  . 
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ure  34.  A  comparison  hot  wen  theory  and  experiment  for  the  log  (  I22)  plotted  against  the  scattering 
le  2  0()  for  the  orientation  Q(  =  90  and  y,=  0  . 
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Figure  35.  A  comparison  betwen  theory  and  experiment  for  the  log 
angle  2  Oq  for  the  onentation  0t  =  90  and  y,=90  . 
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